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Abstract. The energy output in the gamma-ray burst (GRB) prompt emission and afterglow phase is expected to photoionize 
the surrounding medium out to large radii. Cooling of this gas produces line emission, particularly strong in the optical, whose 
variability is a strong diagnostics of the gas density and geometry in the close environment of the burst. We present the results 
of a spectral time series analysis of the host galaxy of GRB 990712 observed up to ~6 years after the burst. We analyze the 
emission line fluxes together with those of the previous observations of the same GRB, in search for photoionization signatures. 
We find that the emission line fluxes show no variation within the uncertainities up to 6 years after the burst, and we use the 
measured line intensities to set a limit on the density of the gas within a few parsecs of the burst location. This is the first time 
that emission from cooling GRB remnants is probed on years time scales. 
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1. Introduction 



In the fireball model of GRBs (e.g. Rees & Meszaros 1992 1, 
the energy released from the collapse of a massive star is con- 
verted into kinetic energy of thin baryonic shells which ex- 
pand at ultra-relativistic speeds. After producing the prompt 
y-ray emission by internal shocks between different shells, the 
residual impacts on the surrounding gas and drives an ultra- 
relativistic shock into the ambient medium. The shock acceler- 
ates relativistic electrons leading to the observed X-ray to radio 
afterglow radiation through synchrotron emission. The emis- 
sion of the X-ray afterglow, integrated over the first 7-10 days, 
typically contains the same energy as the primary y-ray burst 
itself (Rees & Meszaros H9951 . 

Whatever the density of the GRB environment, photoion- 
ization of the circumburst material by the prompt X-ray emis- 
sion and by the X-ray and UV afterglow emission will be in- 
evitable. If the circumburst density is high, it will lead to time 
dependent (on hour timescale) absorption (Perna & Loeb 1998 1 
and emission line features (Bottcher et al. 1999 Ghisellini et al. 
I1999> . such as those claimed to be seen in X-rays. On longer 
timescales, the GRB photoionization may lead to indicative 
recombination line features, which allow the identification of 
remnants of GRBs in nearby galaxies (Band & Hartmann 1992 
Perna, Raymond & Loeb 2000 PRL in the following). 
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Photoionization of the ambient medium is well studied 
in the case of supernovae which are the most similar cases 
to GRBs. For SN 1987A, IUE observations showed that the 
prominent UV lines started to increase simultaneously after 
60-80 days and stayed at a constant level until 400 days after 
the initial exciting supernova outburst (Fransson et al. [T989l. 
After 400 days most lines decreased fastly, and reached the 
noise level by day 1500 (Sonnebom et al. 1997 1. Detailed mod- 
elling of the ionization zones and and the line emission of the 
circumstellar gas (Lundqvist & Fransson 1996 ) allowed to con- 
strain the gas density. 

In the case of a GRB, ionization by the prompt emission, 
the afterglow (photon field) and by the blast wave (shock- 
ionization) will appear additionally to the SN component, and 
largely dominate. The blast wave is expected to influence the 
ionization state of the gas on timescales of hundreds to thou- 
sands of years after the burst (PRL). Therefore, for the pur- 
pose of our analysis, we can safely assume that photoioniza- 
tion is the dominating ionization mechanism. A large fraction 
of the energy (e.g. kinetic, magnetic) stored in the GRB-jet 
(~10 52 erg) is released in the afterglow. The X-ray/UV compo- 
nent of this radiation (with some contribution from the X-ray 
emission of the prompt phase) is the main responsible for the 
ionization of the ambient medium. For typical GRB/afterglow 
luminosities, the size of the ionized region is on the order of 
100 pc for an ISM density of ~ 1 cm -3 (PRL). 

After being ionized, the gas then starts to cool, on a 
timescale t coo \ ~ 10 5 (r/10 5 K)/(n/cirT 3 ) yr, where T is the 
temperature of the gas, and n the electron density. Although 
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cooling starts soon after the burst of radiation has passed 
through the medium, emission from the cooling gas increases 
during the early times after the GRB. If a region of radius R is 
heated and ionized by the burst radiation, the maximum emis- 
sion will occur after a time t ~ R/c ~ 3 yr (R/pc), due to light 
travel times from different parts of the region. 

The strength and timescale of the recombination emission 
depends strongly on the ambient density. While the modelling 
of the broad-band SED of afterglows has led to densities in 
the range 1-10 cm 4 , there are also observational indications 
for much higher densities: (i) observed variable X-ray lines 
(Watson et al. 120021 Reeves et al. 120051 Frontera et al.|2004l 
and continuum absorption (Lazzati & Perna 2002 1 require den- 
sities of ~10 5 -10 6 cm" 3 ; (ii) some GRB afterglow data require 
a dense (~10 4 cm" 3 ) shell around some nearby low-density 
medium (1-10 cm" 3 ) (Chevalier et al. 120031 . The SN-GRB 
connection is now clearly proven for four GRBs (Galama et 
al. [T9981 Hjorfh et al. 120031 Stanek et al. 120051 Malesani et 
al. 120041 Pian et al. [2006 1, indicating the link between long- 
duration GRBs and deaths of massive stars. Observations of 
massive, Wolf-Rayet (WR) like stars, have shown that they 
loose matter via strong stellar winds. A WR stellar wind, in- 
teracting with a circumstellar medium, leads to the formation 
of a shell (termination shock) whose density and radial extent 
depend on both the progenitor characteristics (i.e. mass loss 
rate, wind velocity) as well as the density of the medium (e.g. 
Fryer et al. 120061 . 

An ionized shell of, say, density ~ 10 3 cm" 3 and radial 
extent of a few parsecs, will reach its peak emission on a 
timescale of a few years (due to the light travel times), and 
will cool on a timescale of tens of years. Therefore, the first 
few years after the burst are crucial for detecting cooling emis- 
sion from these dense, compact shells produced by the wind 
termination shocks of the massive stars progenitors of (long) 
GRBs. 

This Letter reports a long-term monitoring, up to 6 years, of 
several strong emission lines, and in particular the OIII/15007 
and the H/J lines, whose ratio is an important discriminant of 
strongly photoionized, cooling gas as expected in GRB rem- 
nants (PRL). Our observations are described in §2, and the re- 
sults are used in §3 to set constraints on the environment of 
GRB 990712. Further implications for studies of galaxy prop- 
erties are discussed in §4, while §5 summarizes our findings. 

2. GRB 990712 

GRB 990712 was discovered by GRBM and WFC onboard 
BeppoSAX on July 12.69655 UT, 1999. The duration of the 
burst was 30s and it was first localized by WFC at R.A. = 
22:31:50, Dec. = -73:24.4 with an error radius of 2' (Heise et 
al. [1999 1. Follow-up observations led to the discovery of the 
GRB afterglow. The redshift of the burst is z=0.433 (Vreeswijk 
et al. 1200 la> . Two different groups found evidence for a SN 
bump from the optical lightcurve of the afterglow (Bjornsson 
et al. I200T1 Zeh et al. 120041 . The BeppoSAX spectra of the 
prompt emission indicated a temporal emission feature located 
around 4.5 keV, which can be fit either with a Gaussian pro- 
file with rest frame energy of around 6.4 keV that is consistent 



Table 1. Log of observations 



Date Instrument Grism Coverage Exposure time 

(nm) (sec) 



14 Jul 1999 


FORS1 


1501 


370 


-770 


2400 


11 Nov 1999 


EFOSC 


Gr6 


400 


-800 


1800 


06 Jun 2002 


FORS1 


600R 


525 


-745 


4320 


17 Jul 2004 


FORS2 


600RI 


512 


-845 


2400 


05-06 Jul 2005 


FORS2 


300V 


445 


-850 


7200 



with an iron line, or with a blackbody spectrum with kT ~ 1.3 
keV (Frontera et al. 2001 1. Since the iron line interpretation re- 
quires a very high density environment that would obscure the 
afterglow, Frontera et al. (2001 1 prefer the thermal component 
interpretation which can be accounted for by the fireball model. 

The host galaxy of GRB 990712 is one of the brightest 
GRB host galaxies with V = 22.3 mag and R = 21.8 mag 
(Sahu et al. 120001 . We obtained spectra of GRB 990712 on 
July 05 and 06, 2005, approximately 6 years after the GRB 
using VLT/FORS2 under good seeing conditions (~0'.'6). The 
300V/GG435 grism/filter was used with a slit width of l'.'O. 
The total exposure time was 2 hours (4x30 minutes). The spec- 
tra were reduced using standard IRAF routines and calibrated 
using standard star G158-100 observed on July 05, 2005 with 
the same grism/filter. The flux calibration was further validated 
by folding the spectrum with the FORS R-band filter curve and 
comparing the obtained magnitude with that of the host galaxy 
given in Christensen et al. (2004). We corrected the spectra for 
foreground extinction of E(B - V) = 0.03 (Schlegel et al. fT998l . 
The line fluxes were determined by fitting a Gaussian to the 
line using the SPLOT task of IRAF. The continuum level was 
determined locally. 

The line fluxes are compared with those derived using VLT 
archival data of observations obtained on July 14, 1999 (PI: 
Galama), November 1 1, 1999 (PI: Courbin), June 06, 2002 (PI: 
Mirabel) and on July 17, 2004 (PI: Le Floc'h) (see Tab[T}. The 
spectra were treated in an identical way, and calibrated using 
standard stars EG 274 (Sep 15, 1999), LTT 377 (at Dec 6, 1999 
and at Jul 15, 2004), and LTT 3854 (Jun 05, 2002) for July 
1999, November 1999, July 2004 and June 2002 data, respec- 
tively. July 1999 data was further flux calibrated using the after- 
glow brightness at that time taken from Sahu et al. (J2000I. The 
flux values we obtain are consistent with the published values 
obtained by flux calibrating the same data using the afterglow 
brightness extrapolated in time (Vreeswijk et al. 12001 a 1. The 
previously unpublished 2002 and 2004 fluxes were similarly 
validated as in the case of 2005 spectra. The estimated error in 
the flux calibration is about 10% for all data. The line fluxes are 
shown in Table [2] The flux errors in Table [2] only include the 
uncertainities in the continuum level. Figure^shows our spec- 
trum from July 2005 and Figure [2] shows the spectrum from 
November 1999 which has the worst signal-to-noise ratio. 

All of the observed lines, which are [O II], [Ne III], Hy, H/? 
and [O III] (/1/14959,5007) lines, are observed to have constant 
fluxes over ~6 years after the burst (see Fig|5}. To overcome the 
difficulty of comparing line fluxes from different spectra ob- 
tained with different settings, i.e. instrument, grism, night con- 
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Table 2. Emission line fluxes 
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Line 






Flux(10- 16 erg/s/cnT 2 ) 






at day 1.5 


at day 123 


at day 1060 


at day 1832 


at day 2185 


[O H](A3727) 


3.52: 


t0.15 




3.66+0.20 




3.40±0.60 


[Nem]U3869) 


0.5: 


tO.l 




0.55+0.07 


0.59±0.05 


0.50±0.05 


Hy(i4340) 


0.3: 


tO.l 




0.53+0.12 


0.46±0.05 


0.46±0.05 


H/3(^4861) 


1.15: 


t0.15 


1.30+0.15 


1.29±0.09 


1.34+0.04 


1.33+0.05 


[0III](/14959) 


2.25: 


t0.15 


2.33+0.18 


2.10±0.15 


2.18±0.08 


2.22±0.08 


[OIII](i5007) 


6.15: 


t0.15 


6.17+0.25 


6.00±0.10 


5.97±0.08 


6.08+0.08 
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Fig. 1. Spectra of the host galaxy of GRB 990712 obtained in 
July 2005. 
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Fig. 2. Spectra of the host galaxy of GRB 990712 obtained in 
November 1999. 



ditions etc., we derive our results on the circumburst environ- 
ment based particularly on the [O III]/15007 and the H/J lines. 
These are close enough in wavelength space to overcome the 
possible effects of different flux calibrations, and at the same 
time, as discussed below, they are especially strong discrimi- 
nants of photoionization vs. collisional ionization. 

3. The circumburst environment 

A distinctive signature of photoionized gas is a high [O III] 
/15007/H/i ratio (k, 5), which is not generally produced in 
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Fig. 3. [Ne III], Hy, U/J, [O II] and [O III] 4959, 5007 line flux 
light curves for GRB 990712. Error bars are smaller than the 
symbol size, if not visible. 



steady shocks for solar abundances (Shull & McKee 1979 1. 
This ratio increases with the ionization parameter and, for typ- 
ical X-ray/UV photoionizing fluxes of GRBs and their after- 
glows, it reaches a value on the order of 100 (PRE). Our obser- 
vations at about 6 years after the burst, as well as the previous 
ones at earlier times, show that [O III] /15007/H/3 ~ 4.6 in the 
host of GRB 990712. Combined with the log ([O II]/H/3) value 
of 0.4+0.1, this value is rather typical of an HII galaxy, as can 
be seen by comparison with the sample of emission line galax- 
ies at redshifts 5= z ^ 0.3 from the Canada-France redshift 
survey (Rola, Terlevich & Terlevich 1997 see also Vreeswijk 
et al. 2001al. Therefore, the constraints that we are able to put 
on the close environment of GRB 990712 will derive from the 
lack of contribution from photoionized gas to the brightest ex- 
pected line ([O III] /15007) from the cooling gas, within the 
observational uncertainties. The measured flux of this line is 
~ 6 x 10~ 16 erg/s/cirT 2 to within ~ 10% uncertainty (see Table 
2). At a redshift of 0.433, this flux corresponds to a luminos- 
ity of ~ 3.8 x 10 41 erg/s for a LCDM cosmology with Hubble 
parameter h=0.73 (Spergel et al. 2006). Since the luminosity 
is constant to within the 10% uncertainty value of ~ 3 x 10 40 
erg/s, the contribution from the cooling GRB remnant cannot 
be larger than this value at any given time during the observa- 
tion window. 

Numerical simulations of cooling GRB remnants (PRE) 
show that the radiation flux from a typical GRB and its after- 
glow ionizes a region on the order of ~ 100 pc for an ISM 
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density of ~ 1 cm 4 . For a gas of solar metallicity, the corre- 
sponding luminosity of the [O III] /15007 line from the cooling 
gas is found to reach a value of ~ 10 38 erg s _1 over a time 
f cooi ~ a f ew xlO 4 yr (see Fig. 3 in PRL; the details of the com- 
putation of the cooling radiation can also be found in PRL). 
The brightness of the line scales with n n e , where n e is the elec- 
tron density. For a gas metallicity not too far from solar (so that 
the particle number is dominated by H) and a highly ionized 
gas, one has n e ~ n, and the line luminosity can be scaled, to 
a first approximation, as 1 ~ 10 38 (n/cirT 3 ) 2 (7? e /100pc) 3 erg 
s _1 . Our observations, up to 6 years after the burst, allow us 
to probe an emitting region, R e , of at most 2 pc in size, due to 
light delay effects as discussed in §1. Since no flux variation 
has been observed within the 10% flux error (corresponding to 
a luminosity of about 3 x 10 40 erg s _1 ), we deduce that the line 
luminosity due to the cooling gas has to be below this level, i.e. 

L5007 ~ 10 38 (n/cnT 3 ) 2 (fl e /100pc) 3 ergs- 1 < 3xl0 40 erg s" 1 (1) 

which yields the limit on the density n S5 6 x 10 3 cirT 3 for 
R e = 2 pc. This limit can be used to constrain the range of 
allowed parameters for the GRB host and ISM densitites. For 
example, Fryer et al. (2006 1 find that, for a progenitor star with 
mass loss M = 1O~ 5 M yr~' and wind velocity ~ 1000 km/s, 
expanding in a medium of density in the range ~ 10 3 — 10 4 
cm 3 , the inner radius of the shell is on the order of tenths of 
a parsec and the outer radius is > 2 pc (with the shell den- 
sity being on the order of the ISM density). These large shell 
densities, filling a region up to the observed emitting volume 
of ~ 2 pc, are not favoured by our observations, since they 
would likely result in a variable [O III] A5Q01 flux over the 
6 years of observations. On the other hand, termination shock 
shells produced in the impact of the wind with a lower den- 
sity medium are consistent with the lack of variability. Also 
note that, given the 6 year timescale of the observations, termi- 
nation shocks located at distances 2 pc cannot be ruled out 
by the currently available data. Generally speaking, observa- 
tions at longer timescales allow one to probe shell terminations 
shocks over a larger range of distances from the GRB progen- 
itor star. However, since more distant termination shocks are 
generally associated with lower ISM densitites (and hence a 
lower luminosity of the [O III] A50Q1 line), in order to separate 
the eventual contribution of this line due to the cooling gas in 
the close GRB environment, from that due to the host galaxy 
itself, a higher signal-to-noise ratio in the observations as well 
as a consistent set of observations (i.e. observations obtained 
with the same instrument and settings) is necessary. 

4. Star formation rate and metallicity 

The star formation rate (SFR) of the host galaxy of GRB 
990712 has been previously calculated based on line fluxes 
and ultra-violet flux, and also an upper limit was derived based 

1 Note however that, as the density increases, the maximum size of 
the region that can be ionized will decrease. Furthermore, note that the 
details of the rise time also depend on the extent of the beaming of the 
ionizing radiation, but we do not worry about secondary effects here 
since the data only allow us to set upper limits. 



on radio non-detection (Hjorth et al. 2000 Vreeswijk et al. 
200 lal EtjOTrJI Christensen et al. I2004i . Using the line fluxes 
of our July 2005 spectrum and the same method used by 
Vreeswijk et al. (1200 la) , we calculated SFR 0// = 2.8+° 4 M Q /yr 
and Ay = 1.7+^ mag based on the Hy/HyS ratio. Our Ay value 
is on average lower than the Ay = 3.4+ 2 4 derived by Vreeswijk 
et al. (200 la 1, therefore our extinction corrected SFR value 
{SFRoihorr = 10+^ 5 Mg/yr) is also lower than their calculation. 
Nonetheless both Ay and the extinction corrected SFR are in 
agreement with Vreeswijk et al. (2001a 2001b 1 values within 
the errors. 

We do not detect the [O III] /14363 line which is neces- 
sary for the determination of the electron temperature T c by 
means of lines (like [O III], [Ne III] etc.) from high ionization 
zone elements. Therefore, we cannot estimate the oxygen and 
neon abundance based on the electron temperature. However 
we can still estimate the oxygen abundance based on the ra- 
tio R 32 = ([OIIM3727 + [0111^4959,5007)/^. Kewley & 
Dopita (2002 1 suggest to use the relation given by Zaritsky 
et al. (1994 1 to obtain an estimate of the oxygen abundance. 
Using the formulae given in Zaritsky et al. ( 1994 1 we obtain 
log(0/H)= -3.7+0.1. However, that formula is calibrated for 
metal rich galaxies and overestimates the metallicity for val- 
ues log(0/H) < -3.5 (see Kobulnicky et al. IT9991 Kewley & 
Dopita 2002 1. Therefore we used equation (16) of Kobulnicky 
& Kewley (2004i, which is adopted from the relation given 
by Kewley & Dopita (2002 1 and parameterized for the lower 
metallicity branch (log(0/H) < -3.6). The result is log(0/H)= 
-3.7+0.1, which is the same as our initial estimate using the re- 
lation given by Zaritsky et al. ( 1994 1. Similarly, Vreeswijk et al. 
(2001 a 1 obtained -3.7+0.4 for log (O/H), which is in agreement 
with our estimate. The oxygen abundance we obtained is just a 
bit lower than the solar value (log(0/H) = -3.34; e.g. Asplund 
et al. 2005). Therefore our assumption that the metallicity of 
the gas nearby the GRB is not far from solar is reasonable in 
our derivation of the [O III] line luminosity of the cooling gas. 

It should be noted how, in the context of GRBs, studying 
the photoionization signatures of cooling gas not only helps 
understand the nature of the GRB progenitor star, but it also 
helps reduce possible biases in the determinations of two im- 
portant properties of the GRB host galaxy such as the SFR and 
the metallicity. In fact, SFR and metallicity calculations based 
on emission lines generally rely on observations obtained at a 
single epoch. However, in order to properly assess the possible 
level of contamination by the cooling GRB remnant, multiple 
epochs of observations spanned over a long timescale are nec- 
essary. 

5. Summary 

We have presented the results of the spectral analysis of the 
host galaxy of GRB 990712. With the last set of observations 
taken about 6 years after the burst, this is the longest time cov- 
erage for a GRB host galaxy up to date. Though we do not 
detect line variations, timescales of a few years are important 
for detecting cooling radiation from the heated shells produced 
by the wind termination shocks of the massive star progenitors 
of the GRBs. For the case of the GRB 990712 host, the lack 
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of time variability in the [O III] /15007 line, combined with the 
^ 5 ratio of the [O III]/l5007/Hj8 lines, has allowed us to set an 
upper limit to the contribution from the cooling gas. This limit, 
in turn, could be used to constrain the allowed range of densi- 
tites within a region of about 2 pc surrounding the burst. We 
have therefore shown how this type of observations provides a 
useful complement to the studies of the close environment of 
GRB progenitors and, therefore, can help reconstruct the char- 
acteristics of the GRB progenitor star. 

Finally, we have pointed out how, if a substantial contam- 
ination to the galaxy spectra is provided by the GRB cooling 
radiation, inferences of the SFR and metallicity that are drawn 
from measurements of line ratios can be biased. To be able to 
assess the degree of this contamination, long-term monitoring 
of the GRB host galaxies is necessary. 
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